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Abstract 

Two photon collisions at high energy have an important theoretical advantage: the 
simplicity of the initial state, which gives us a unique opportunity to calculate these pro¬ 
cesses for large virtualities of both photons in perturbative QCD approach. In this paper 
we study QCD saturation in two photon collisions in the framework of the Glauber-Mueller 
approach. The Glauber-Mueller formula is derived emphasizing the impact parameter de¬ 
pendence (bt) of the dipole-dipole amplitude. It is shown that non-perturbative QCD 
contributions are needed to describe large 6t-behaviour, and the way how to deal with 
them is suggested. Our approach can be viewed as the model for the saturation in which 
the entire impact parameter dependence is determined by the initial conditions. The uni- 
tarity bound for the total cross section, its energy dependence as well as predictions for 
future experiments are discussed. 

It is argued that the total cross section increases faster than any power of ln(l/x) in a 
wide range of energy or x , namely <7(7* —7*) oc (1/Q^) exp{a^/ln{l/x)) < l/mf where 
exp{ayjhi{l/x)) reflects the x dependence of the gluon density xG oc exp{2a^/\nf\Jx)) 
and tUtt is the pion mass. 
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1 Introduction. 


Two photon collisions at high energy have three theoretical advantages over hadronic collisions 
or/and deep inelastic scattering: 


• The simplicity of the initial state, which allows processes, such as large transverse mo¬ 
mentum hadronic jet production, to be calculated exactly to lowest order in perturbative 
theory. With the advent of high quality experimental data, theoretical analysis also focus 
on higher order corrections to the basic processes which can provide an interesting test 
of the theoryP; 

• Scattering of two photons with large but equal virtualities gives unique access to BFKL 
emission j2j, making this process very useful for studying the dynamics p in 13 mu El i; 

• Scattering of two virtual photons with large virtualities allows one to study shadowing 
(screening) corrections on the solid theoretical basis of perturbative QCD |l()j . 


It is well known that the correct degrees of freedom at high energy are not quarks or gluon 
but colour dipoles nn ca ESI cn which have transverse sizes and the fraction of energy 
X. Therefore, two photon interactions occur in two successive steps. First, each of virtual 
photon decays into a colour dipole ( quark - antiquark pair ) with size r*. At large value of 
virtualities the probability of such a decay can be calculated in pQCD. The second stage is the 
interaction of colour dipoles with each other. The simple formula ( see for example Ref. ini 
) that describes the process of interaction of two photons with virtualities Qi and Q 2 is (see 

Fig.E) 

Nf 

dXY. (I'l) 

a,b 


a{QuQ2,W) = I 


dzi J ^ dz2 j d^dr2,t|T^_^(Q2; 2 ^ 2 , ^2,01^ o-aJ(a^ab, n,*, ra,*; 6*) 

where the indices a and b specify the flavours of interacting quarks, T and L indicate the 
polarization of the interacting photons. The r* denote the transverse separation between quark 
and antiquark in the dipole ( dipole size) and Zi are the energy fractions of the quark in the 


photon i. 


cr 


dd 

ab 


= 2N{{xab,ri^t,r2^t',bt) where N is the imaginary part of the dipole - dipole 


amplitude at energy x given by 


~ _ Ql + Ql + + Ami 

+ Qj + Ql ^ ^ ^ 

where is the mass of the quark with flavour a. bt is the impact parameter for dipole-dipole 
interaction and it is equal the transverse distance between the dipole centers of mass. It is clear 
that afb has a meaning of da/d?bf 

The wave functions for virtual photon is known Hi and they are given by 
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Zl + {I- zy)QlK^{Qart) + mlKl{Qart)) (1.3) 


(1.4) 


with = z{l — z)Q‘^ + where Za and nia denote the faction of charge and mass of the 
quark of flavour a. 



Figure 1: The picture of interac¬ 
tion of two photons with virtu¬ 
alities Qi and Q 2 larger than a 


’soft” scale. 


V(r,Q') 


The main contribution in Eq. (HH) is concentrated at ri^t ~ 1/Qi ^ 1/h ri^t ~ 
I/Q 2 -C l//i where fi is the soft mass scale. Therefore, at hrst sight, we can safely use pQCD 
for calculation of the dipole-dipole cross section a in Eq. (HH). The objective of this paper is 
to investigate the dipole-dipole cross section at high energy ( low x ) where QCD saturation 
is expected naiiHiiini. The hrst analysis based on Golec-Biernat and Wiisthoff model [20] 
has been performed in Ref. uni. Here we will extend this analysis by using the Glauber - 
Mueller approach uniniiioi with special focus on the impact parameter dependence which 
was completely omitted in the GBW model as well as in Ref. HQ. 

In the next section we discuss the dipole-dipole interaction in the Born approximation of 
pQGD. We show that this approximation leads to a which decreases as a power of bt- It turns 
out that a 1/6^ for large bt > ri^t and r 2 ,t. Of course, such a behaviour will not change its 
character in higher orders of pQGD ( see Refs I2III221I23E3 ) since it is a direct consequence of 
the massless gluon in QGD. Using the Born approximation as the example we consider the non- 
perturbative contribution that provides an exponential decrease at large values of bt > l/rriT^. 

Section 3 is devoted to Glauber - Mueller formula in the case of the DGLAP emission 
I2ni- Here, we use the advantage of photon - photon scattering with large photon virtualities, 
since we can calculate the gluon density without uncertainties related to non-perturbative 
initial distributions in hadronic target. It is well known that no bt dependence is induced by 
DGLAP emission at least for large values of the impact parameter. Therefore, the entire impact 
parameter dependence is due to the Born Approximation cross section. In other words, we can 
use our approach as an explicit illustration of the point of view that the non-perturbative large 
bt > where is the pion mass, is determined by the initial condition (see Refs. 
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Figure 2: The Born approximation for dip ole-dip ole scattering amplitude. 


121 ED]) in the contrast with the notion that such a behaviour could change the kernel of the 
non-linear equation that governs evolution in the saturation region 12Ill2ai2SlE3- 

The unitarity bounds as well as different regimes of the energy behaviour of two photon 
total cross sections are considered in section 4. 

In section 5 we give our estimates for the values of the total cross sections for accessible 
range of energies. 

In the last section we summarize our results. 


2 Dipole-dipole interaction in the Born approximation. 


The Born approximation for the dipole-dipole scattering amplitude is shown in Fig. |21 

To obtain the expression for ri r 2 ^tt] h) (see Fig.[T]) we need to calculate the diagrams 

in the momentum representation and than to rewrite them in space-time representations. The 
conjugated variables to pt and k will be the size of the dipole (say ri^t and the impact parameter 
bf. The detailed calculation performed in light-cone technique ( see for example Ref. [211 ) has 
been performed in 123. The answer is 




^ 2 Nl-I ( - Zin - Z2r2y ip - Zin - Z2r2Y \ 

2iV2 ^ (^5 _ _ Z2r2y {b - zifi - Z2r2Y J 


(2.5) 


where Zi is the fraction of the energy of the dipole carried by quarks and Zi = Zi — 1. All 
vectors are two dimensional in Eq. (j23|). 
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Eq. (j2.5|l has a simpler form if we assume that Zi = 1/2. Namely, 


o{x,rx^ur2,tM) 


7ral 


2 AT/ 


(^^^ jb + RY {b-Ry \ 
[ (b + S)2 (6 - S)2 ) 


( 2 . 6 ) 


where R = ^ and S = We note that we do not hnd the dipole-dipole cross 

section in impact parameter representation in Ref. m, but the calculation is so simple that we 
just present the answer. 

To simplify our further calculations we restrict ourselves by DGLAP emission assuming that 
ri t is much smaller than r 2 ^f If is instructive to hnd two different limits in Eq. (ESD- 


• > 

angle 


f' 2 ,t S> Expanding Eq. (Eini) one can obtain after integration over azimuthal 


(T{x,ri^t,r2,t]bt) 


2^c -1 
c 


' l,t ' 2,t 

bf 


(2.7) 


• bt ^ ri^t < r2,t- 


(T{x,ri^t,r2Rbt) 


ira 


2 N!-i 

iV? 


r‘2 r‘2 

' l,t ' 2,t 

^2 

^■2^2 ' 2,t 


( 2 . 8 ) 


• Therefore, we can suggest a simple formula which covers two these limits, namely 


(y{x,rx,ur2,ubt) 




2N^c 


' 2,t 


«? [ hi + Z 2 Z 2 r|, )= 


(2.9) 


For further estimates we will use Eq. (El) which rehects all qualitative features of the 
full expression of Eq. (EH but considerably simplifies the calculations. 


Eq. fj2.5p as well as Eq. (EH leads to power - like decrease at large values of bt, namely, 
cr(5;, ri^t, r 2 ,p 6i) oc Such behaviour cannot be correct since it contradicts to the general 
consequence of analyticity and crossing symmetry of the scattering amplitude. Since the spec¬ 
trum of hadrons has no particles with mass zero, the scattering amplitude should decrease as 
g- 2 m,r bt _ Certainly we need to take into account non-perturbative corrections to heal this 
problem as has been noticed in many papers [21 ED El El El ED]. We suggest the procedure 
how to introduce such corrections which is based on the hadron-parton duality in the spirit of 
the QCD sum rules m. This procedure consists of two steps: (i) hrst, we rewrite Eq. EH) 
in the momentum transfer representation [t = —ql) in the form of a dispersion relation with 
respect to the mass of two gluons in f-channel; (ii) secondly, we claim that this dispersion 
integral gives correct contribution of all hadronic states on average. Therefore, the model for 
the non-perturbative contribution is the integral over a two gluon state in the f-channel but 
with the restriction that two gluon mass should be larger than the minimum mass in hadronic 
states, namely, larger than 2m^. As in QCD sum rules EH we assume that the integrand at 
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large mass of two gluon state can be found in perturbative QCD, while for small mass we have 
to include the realistic (experimental ) spectrum of hadrons. The integration from 2171^^ means 
that we believe that we can approximate the dispersion integral even in region of small masses 
by the perturbative QCD contribution. This procedure corresponds to the approximation that 
has been used in Ref. li]|. We can also evaluate this integral differently: to take into account 
the hrst resonance ( glueball ) explicitly and to use pQCD approach to estimate the dispersion 
integral for masses larger than sq, the value for sq can be taken from a QCD sum rules calcula¬ 
tion of the glueball spectrum [21]. Such an approach is closely related to one developed in Ref. 
1221 and it appears reasonable in pure gluodynamics where we do not have any pions. 

Rewriting Eq. (I23D in the form 


(T{x,ri^t,r2,t;bt) 


with obvious notation, we can see that 


(62 + a2)2 ’ 


= CK„r,,) / b, 


tdbt 


Mbq) 
(6? + 


= C(rij,r2f) ^Ki(aq) 


( 2 . 10 ) 


( 2 . 11 ) 


where Jq and Ki are Bessel and McDonald functions respectively. However, we can rewrite 
Ki (a q) in a different way as 


qKi{aq) = j 


Ji(k a) dn 


+ q^ 

The last integral (see Eq. (Em) gives as the dispersion relation, namely, 

1 Ji{Ka) IT? dK 


1 

(^{x,ri^t,r2,t]t =-q^) = C(ri,i, rs,*) — 

2a Jo 

Eq. ()2.1djl we replace by 


H? — t 


= -q^) = C{ri^t,r2,t) — [ 

2a J 2 


1 /■°° Ji{Ka)K?dK 


q t? - t 


( 2 . 12 ) 


(2.13) 


(2.14) 


accordingly to our main idea. Returning to the impact parameter representation we obtain 


o'(a:,ri,t,r2,t;6t) 


poo poo 

<^(^,*,^ 2 ,*) — / n'^dn qdq 
2a 22m,r Jo 

2 - 1 2 2 1 r 

^ 2iV2 ’ ’ y/Z2 Z2 72,t J2i 


Ji{Ka) Jpjqbt) 

XD 

dn Ji^Ko) KQ^Kbt) .(2.15) 


One can see that a oc e for bt S> 1/2 due to the asymptotic behaviour of the 

McDonald function Ki{Kbt) —> at large bt- 

Therefore, the bt behaviour is: for 1/(2 m^) > bt > ri ^ or/and r 2 ^t the dipole-dipole 
scattering amplitude falls as l/bf, but for large bt {bt > l/{2m.,^)) we have normal exponential 


6 













(DGLAP emission) (Giauber - Mueiier formuia) 


Figure 3: The Glauber - Mueller approach for dip ole-dip ole scattering amplitude. 

decrease as which has a non-perturbative origin. Eq. (ITT31) gives us a rather general 

way to take into account the non-perturbative contribution, since in this equation we explicitly 
introduce the minimum mass in the experimental hadronic spectrum. However, as we have 
mentioned above, we can expect a large mass for the low limit of integration in the dispersion 
relation of Eq. (imi) cLnd Eq. fj2.15j) (Qo ^ ^tt) which will lc8jci to bchcLviour 

oc 


3 Glauber - Mueller formula. 

Glauber - Mueller approach takes into account the interaction of many parton showers with 
the target as it is shown in Fig. 121 Actually this formula was suggested in Refs. [I21II21 but 
Mueller im was the hrst who proved this formula especially for gluon parton density. The 
main idea of this approach is that colour dipole is the correct degree of freedom for high energy 
scattering h Indeed, the change of the value of the dipole size r* (Ar^) during the passage of 
the colour dipole through the target is proportional to the number of rescatterings (or the size 
of the target R) multiplied by the angle kt/E where E is the energy of the dipole and kt is the 
transverse momentum of the t-channel gluon which is emitted by the fast dipole. 

Ar^ oc ■ (3.16) 

Since kt and rt are conjugate variables and due to the uncertainty principle 

kt (X — . 

rt 

Therefore, 

k. 1 

Art oc R— <C rt if R rf E or x -C -=: . (3-17) 

E 2m R 

Mhis idea was formulated by A.H. Mueller in Ref. m a bit later. 
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3.1 DGLAP emission 


We first discuss the generalization of the Born approximation to include the DGLAP emission 
( see Fig. IHl). This will give us the correct description of the one parton shower interaction. 
The DGLAP equation looks very simple in the region of low x, namely, 


d^xG{x,rlt,rlt) 
dhi{l/x) d ln(l/ri j) 


xG{x,rl^, 

TT 



(3.18) 


where we consider ri^t "C r 2 ,t and rewrite the DGLAP equation in coordinate space. We 
would like to recall that the DGLAP evolution equation sums the (05 log contribution, 
and therefore, we can safely rewrite it in the coordinate representation since within logarithmic 
accuracy InQ^ = ln(l/r^). The initial condition for Eq. (jd.lHj) is xG(x = xo,rf^,r 2 t) = 1. It 
means that the dipole-dipole cross section at hxed bt for one parton shower interaction has a 
form ( 

(rdipoie{x,ri^t,r2,t]bt) = xG{x,rlt,rl^) , (3.19) 

where (Tfi^oie is the Born approximation for the dipole cross section. 

The obvious solution is 

xG{x,rlt,rl^) = Jo (^ 2 ^^(ri,t,r 2 ,t) ln(l/a;)^ (3.20) 


where 


^(D,t,r2,t) 


12 AT, ln(4/(r2^A2)) 

11W - 2Nf ''ln(4/(r2%A2)) ' 


Here, in the arguments of the running QGD coupling we have made the simple replacement 
—>■ 4/rj. Within log accuracy we cannot guarantee the coefficient 4 in this expression but 
as was argued in Ref. [SI] this is a reasonable choice (approximation). 

Eq. (13.2011 has the following asymptotic behaviour 


xG{x,rl^,rl^) 


g2 y/Ciri,t,r 2 ,t) \n{l/x) 


(3.21) 


which means that xG grows faster than any power of ln(l/x). 

Strictly speaking Eq. (ITTHD is proven in so called double log approximation of perturbative 
QGD, in which we consider 


as^^{^lx)\n{Ql/Ql) ^ 1 ; 

Q!5ln(l/a;) < 1; 

as\n{Ql/Ql) < 1; 

as -C 1. 


(3.22) 


However, we will use this equation in a wider kinematic region where 0^5 ln(l/a;) \a{Q\/ Q^) > 1 
and Q( 5 ln(l/a;) > 1 while as^^iQilQl) ~ kinematic region we should use the 
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BFKL equation [2]. We view Eq. (IXTHl) as the limit of the BFKL equation in which we take 
into account the logarithmic contribution in transverse momentum integration in the BFKL 
kernel. The justihcation for such an approach is the fact that the anomalous dimension of the 
DGLAP equation can be parameterized in simple way [22l 

7M = - l) (3.23) 

The hrst term in Eq. leads to Eq. ()3.18|) in the region of low x. 

Eq. (I3.1hj) and Eq. (I3.2()|l solve the problem of one parton shower interaction in the DGLAP 
evolution. It should be stressed that the impact parameter dependence enters only in the Born 
term in Eq. (inni). The explanation of this fact is very simple if we recall that the logarithmic 
contribution originates from the integration over transverse momenta large g^/4 where is 
the momentum transfered along the DGLAP ladder. Therefore, we have two choices: (i) 
t > Q 2 and in this case all logs can be summed in a function with the argument ln(Q^/g^), (ii) 
t < Q 2 and in this case we have function of as in Eq. (I3.2()j) . In our problem we are 

certainly dealing with the second case since we are mostly interested in large bt - behaviour of 
the scattering amplitude which corresponds to low behaviour. 


3.2 Many parton showers interactions 


Since the colour dipoles are correct degrees of freedom the unitarity constraints are for dipole- 
dipole elastic amplitude aei{x,ri^t,r 2 ,t',bt) are diagonal and they have the form 


2 Imaei{x,ri^t,r2,ubt) = (3.24) 

a{x,ri^t,r2,t;bt) = |aei(^) ^2,d ^t)P + Gin{x,ri^t,r2,t;bt) 

where Gin stands for contribution of all inelastic processes. Eq. is exact for dip ole-dip ole 

scattering while it has only limited accuracy, for example, for dipole-proton scattering (see Ref. 
m )• An experimental manifestation of the poor accuracy of Eq. for deep inelastic 

scattering is the large cross section of so called inelastic diffraction dissociation of proton in an 
excited state. 

Assuming that at high energies the amplitude is pure imaginary, one can End a simple 
solution to Eq. (IT^ . namely, 

aei{x,ri^t,r2,t;bt) = i (1 - e 2 \ ■ (3.25) 

GU^,n,t,r2,t;bi) = (1 - ; (3.26) 

where D is the arbitrary real function. 

In Glauber - Mueller approach the opacity D is chosen as 


^{x,ri^t,r2,ubt) 


(^di^oie{x,ri,t,r2,ubt) 


(T. 


diioie{x,ri,t,r2,t-,bt) xG{x,rlt,rlt) , 


(3.27) 
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where is dipole-dipole cross section in the one parton shower approximation (see Fig. |2)). 

One can guess that the physical interpretation of Glauber-Mueller formula is simple, namely, 
it takes into account many parton shower interactions in dipole-dipole scattering, but it does 
not include a possibility for produced partons from different parton showers to interact. These 
interactions lead to a more complicated non-linear evolution equation (see Refs. [miiHiiiniEH] 
). The influence of non-linear evolution on the photon-photon scattering will be discussed in a 
separate publication; here we restrict ourselves to consider only the hrst step of this non-linear 
evolution which is the Glauber-Mueller approach. 


4 Unitarity bound 


Using the Glauber - Mueller formula of Eq. (FOKli we can give the unitarity bound for dipole- 
dipole scattering as well as for 7 * — 7 * total cross section (see Eq. (HH)). We consider the 
Glauber - Mueller formula for the total dipole-dipole cross section, namely. 


a 


dd 

tot 



2 


(4.28) 


with opacity G is given in Eq. (HOTI) . 

The main idea m is to replace the full integration over impact parameter in the expression 
for the total cross section, by integration in two different regions: (i) the hrst region is 0 < 
bt < bQ{x)] and (ii) the second one bo < bt < cx). In the hrst region we consider Q/2 > 1 
and replace Im aei by 1. On the other hand, in the second region we assume that 0/2 < 1 and 
expand Eq. with respect to O restricting ourselves to the hrst term of this expansion. 

Therefore, 

ffS < 2ir (^//“ dbl 1 + /r l) ■ 



Let us assume that bo -C 1/2 . In this case we can use Eq. ()2.9|1 (or Eq. ()2.10|l l for the 

bt - dependence for both intervals. Taking the integral of Eq. we have 


a 


dd 

tot 


< 2n 


boix) 


+ 


C{ri^t,r2,t) xG{x,rl^,rl^) 
2{bl{x)+a^) 


(4.30) 


We follow Froissart’s idea to evaluate the value of the characteristic impact parameter bo, 
namely, the value of bo{x) can be found from the following equation 


^{x,ri^t,r2,ubo{x)) 

2 


(4.31) 
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Indeed, for bt > bo ^ < 1 the full formula gives less than the hrst term of the expansion, while 
for bt < bo ^ > 1 , the elastic amplitude for the hxed value of the impact parameter is less 
than 1. Using Eq. (Id.lDj) we obtain the solution of Eq. (I4.dlj) in the form 


bl = -rlt{z2Z2) + ln(l/3;)^ (4.32) 

see Eq. for all notation. 

One can see that Eq. (I4.32j] leads to b^ — >• in{i/a:) ^ g means that 

bo{x) increases faster than any power of ln(l/x). 

One can see that 6 q becomes negative at rather large values of x. It reflects the fact that 
Eq. (jd.dlj) does not have a solution at all values of x. In other words, 0/2 < 1 even at 6 * = 0 
for low energies. However, it should be stressed that Eq. (lOTl) does have a solution at high 
energies which we are actually dealing with in this paper. 

Substituting Eq. (lO^ into Eq. (iraii we obtain 


< 271 (^2bl{x) + rl,(z2Z2)) 


(4.33) 


with bo of Eq. (lO^ . Eq. (lO^ is in striking contradiction with the Froissart theorem which 
states that ln^(l/x) (see R ef. for more details on the Froissart theorem for photon 

interaction.) 


bo » 


2 m,r 


In this case for bt > bo > 2 ^ ^^e integral over k in Eq. (ITT^ is concentrated at k —>■ 2 
with K —2m^ ~ 1/6*. Since 2m^a 1 we expand Ji function, namely Ji^na) = ^/ta. 
Since we are interested only in large bt behaviour of the opacity H (60 3> y^) we replace 
under integral by (2m^)^. The use of the asymptotic behaviour of McDonald’s function as well 
as the simplihcations, mentioned above, leads to overall accuracy 1/bt in the pre-exponential 
factor, which is enough to obtain the unitarity bound. It worthwhile mentioning that in our 
numerical calculation the integral of Eq. was computed without any approximation. 


Finally we have the following estimate for the integral of Eq. dnni)^ 


BA 

^ dipole 


2 - 1 

2iV/ 

o N^-l 




POO 

'm ^ 2 ,t )(2"1^)M ndn Koinbt) ] 

’ ’ J2 m-rz 

(4.34) 

M rlt){2m^f^ Ki{2m^ h) ; 
bt 

(4.35) 


(4.36) 


§The integral zKo{z)dz = Ki{2nn.,r). It follows directly from the differential equation for Kq, namely 
-^{zKo{z)) = —zKo{z) which should be integrated over z. Recalling that —-^Ko{z) = Ki{z) we obtain the 
above integral. The alternative way is to use a combination of equations 6.561(8) and 6.561(16) from Ref. |d9| . 
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where Eq. (j4.H(i|l gives the asymptotic behaviour at large bt {pt ^ l/{2m.P) ). Namely, this 
is the expression we will use for the estimates of the value of 6o in this case. Substituting 
Eq. (l4.Hb|l in Eq. (I2.15j) and Eq. we hud the solution to Eq. (lOTl) . namely 




j J '■xG{x,rlt,rlp] (4.37) 

\/4m^6o’ ^{x) J 


Eq. ()4.37|1 is still an equation for 6“^ which has the asymptotic solution at low x: 


= 


1 , ( 2 - 1 , 2 2 
In { TV a ( 


2 


Am 




(4.38) 


One can see two important differences between this case and the case that we have considered 
previously : 

1. &o of Eq. (I4.38|l grows only logarithmically as a function of energy. From Eq. (I3.2()j) we 
conclude that oc Y^ln(l/a;); 

2. the second term in Eq. gives a small contribution which does not depend on 

energy. 


Therefore, in this kinematic region the unitarity bound has a form 

^dipole-dipole ^ 27l{bl^P{x)f ( 4 . 39 ) 


with 6“^ from Eq. (j4.38jl . 

This equation reproduces the classical Froissart result [2H], namely, the fact that the total 
cross section can increase only logarithmically. This is the kind of energy behaviour we expect 
for DIS or hadron-hadron collisions. However, we would like to draw your attention to the fact 
that we obtain 

Ptt 

Hl/x) ; ( 4 . 40 ) 

while the unitarity bound for the hadron - hadron cross section has In^ s - behaviour < 


It is worthwhile mentioning that Eq. (lonii holds in the wide range of the photon virtualities 
which we will dehne below. 

Predictions: 

Comparing Eq. (lOHll and Eq. (lonn . one can see that in a wide range of energies where 

boi^x) < 1/2 tUtt the photon - photon scattering shows an exponential (oc ^ behaviour 

as a function of ln(l/a;), in striking contradiction with the DIS or/and hadronic processes. 
However, for higher energies b^Q^{x) reaches the value of l/2m7r or 1/mgiy^ebaii- For higher 
energies the unitarity bound becomes the one of Eq. (irail . The numerical evaluation shown in 


12 



























y=log,(|(1/x) 


Figure 4: Energy behaviour of bo{x) ( = 4:/r‘f). 

Fig. |3] illustrate the fact that the kinematic region of an exponential increase is wide, especially 
if we believe that the non-perturbative corrections will only appear at small masses in t-channel. 
Therefore, we hnd that 7 *— 7 * scattering shows quite different behaviour than DIS and hadronic 
processes at all accessible energies (see Fig. ID). However, if the typical mass in t-channel is 
rather the mass of a glueball (see Ref. jl2l ) fhe non-perturbative corrections will stop the 

exponential increase as at x ~ 10 “^. 

It is interesting to notice that the value of bo{x) turns out to be larger at larger value of 
in the region of low x. The reason for such behaviour is the fast increase of the gluon density at 
larger values of which prevails the suppression due to extra factor 1/Q in Eq. (I4.d2j) . From 
Fig. Hlone can see that &o(Q^ = 20GeV‘^) < 6 o(Q^ = 40GeR^) at x < 10“^. 
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5 Total 7 * — 7 * for accessible energies 


Using the master formula of Eq. (HID with the dipole-dipole cross section given by Eq. 
we calculate the 7 * — 7 * total cross sections at accessible range of energies. The results of 
calculations are presented in Fig. El We £x the virtuality of one of the photons at Q 2 = 4GeU^ 
and calculate the cross section at different values of Q^. It is essential to recall that we discuss 
7 * — 7 * scatterinbg in the DGLAP dynamics and we have to £x large values for virtualities 
of both photons. Q 2 = 4GeU^ corresponds to r 2 ^t ~ 0.2 fm which is smaller than the 
electromagnetic radius of pion = 0.66 fm). Therefore, we can apply perturbative QCD to 
our process. 

In Fig. 0 we also show the experimental data for the 7 * — 7 -process since there is no 
experimental information about the values of the cross sections for 7 * — 7 * - scattering for large 
but different photon virtualities. However, the main dependence of the cross section is on the 
largest virtualities and we can hope that the data on 7 * — 7 reaction is not very different from 
7 * — 7 * one. 

One can see from Fig.0that our predictions are not in contradiction with available but poor 
experimental data. We see in Fig. 0 that data with one real photon overshoot our predictions. 
Actually, there are more data on 7 * — 7 reaction but they are presented in the form of photon 
structure function. We do not want to recalculate the cross section using these data since we, 
being theorists, are not entitled to put experimental errors for these reconstructed data. We 
would like to mention once more that this comparison with the experiment could be considered 
only as illustrative one showing that we obtain a reasonable estimates for the value of the cross 
sections. The fact that the data with large but equal virtualities are less than our prediction 
is understandable since in our approach the cross sections for such processes do not have an 
extra enhancement due to the gluon structure function. 

Therefore, we can view Fig. 0 as an argument that our predictions do not contradict the 
current experimental data and as the reason for our expectations that future experiments will 
provide us with data which we will be able to compare with our predictions. 

It should be mentioned that for serious comparison with the experimental data we have to 
calculate the power-like corrections to high energy behaviour discussed in this paper. These 
corrections are calculable for 7 * — 7 * scattering and they can be described as the exchange of 
quark-antiquark pair in f-channel ( so called “box” diagram of Fig. ED EH. The simple “box” 
diagram without gluon emission falls down as 1/lU^ where W is the energy of 7 *— 7 * scattering. 
However, the gluon emission slow down this decrease and, therefore, such corrections could be 
important at sufficiently high energies. 

It is instructive also to compare the realistic calculation with the unitarity bound (see 
Fig. ED. To calculate the unitarity bound we use Eq. (ED where we substitute 

/ = 2ti (^2hl{x) + rl^{z2Z2)) . (5.41) 

One can see that the unitarity bound considerably overestimates the value of the cross 
section. 
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Figure 6: The picture of in¬ 
teraction of two photons with 
virtualities Qi and Q 2 due to 
quark-antiquark pair exchange ( 
so called “box” diagram). 
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Figure 7: Energy behaviour of total 7 * — 7 * cross section for low energies and experimental 
data. Squares denote the L3 data [33] while the triangles mark OPAL data [331 • Circles label 
data taken from E|. 
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Figure 8 : Energy behaviour of total 7 * — 7 * cross section (solid line) and unitarity bound 
(dotted line). 

6 Summary 

We can summarize our approach in the following way. The kinematic region which we study 
in this paper is the high density QCD region. In this region we have the system of partons at 
short distances at which as is small, but the density of partons has become so large that we 
can not apply the usual methods of pQCD. The important method to deal with hdQCD is a 
Glauber- Mueller approach, which gives the simplest approximation for the high parton density 
effects. Developing the Glauber- Mueller approach, we obtained the following results. 

• Both DGLAP and BFKL equations are linear evolution equations predicting a steep 
growth of cross sections as a function of energy. However, it is believed that unitarity 
holds for all physical processes. At high energies it manifests itself as a suppression of the 
growth of the cross section. At the saturation scale Qs{x) nonlinear effects set in. These 
effects are due to formation of a high density parton system. 

• In this paper for the hrst time the Glauber - Mueller approach has been developed for 
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the case of virtual photon - photon scattering. It allows us to estimate the saturation 
scale where the transition occurs from the low density to the high parton density regime 
. The estimate is made from the equation ( ri > r 2 ) 


^{b ^ 1'f'l,saturation) '^ 2 ) 1 (6.42) 

The solution of this equation is shown in Fig. El 



\og^^ (1/x) 


Figure 9: Estimate of the saturation radius where the transition occurs from low density to 
high parton density regime. Solution given as dependence of r'^aturation fixed 

values of Q^- 


The solution to Eq. 0iL42|) is proportional to rl^atuvation oc (a:G(x, j) \ 

It is not surprising that the value of ri^saturation decreases as function of Q 2 as one can 
see in Fig. El At first sight it looks strange that the value of the saturation scale Ql = 
’’^‘saturation rather large. Indeed, for x = 10^ and Q 2 = 4Geld^ the value of Ql ~ 
70 GeV^ is much larger than expected ~ 1 — 2 GeV^ for proton. 
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To understand this difference we take the parameterization of the gluon structure function 
for proton in the form of Eq. (Id. 2 ()|l .namely xG{x, = GoIo{2 Rp) ln(l/a;)) [Hlj . 

Rp is a proton radius which in this estimate we can take R^ = lQGeV~^. Go is equal 
0.136. One can obtain 


~ proton) 


cx QlRlil/Go) ^ 70. 


Therefore, we claim that the large value of the saturation momentum is one of the inter¬ 
esting features of the 7 * — 7 * scattering at high energy. 

• We note that the gluon interaction leads to power-like decrease of the opacity (O) in 

Glauber - Mueller formula as a function of the impact parameter (b), namely G oc l/bf. 
It turns out that because of this behavior the 7 * — 7 * cross section has a wide range of 
energy where it increases faster than any power of ln(l/x) in remarkable contradiction 
with hadron - hadron and deep inelastic cross sections, which cross sections can have only 
In^ W growth with energy [2H] . This fast rise of the 7 * — 7 * cross section continue up 
to energies at which the typical impact parameter ( 60 ( 2 ^)) will reach the value of 1 / 2 mjr 
{bo = 1/2 , see Fig. Hj). 

• The influence of this power-like bt behaviour on unitarity bound is studied. This bound 
is calculated to give an estimate for the energy behavior of the cross section. 

• It is shown that non-perturbative contributions are needed even for the case of photon- 
photon scattering with large virtualities of both photons in order to describe the large 
6 t-behavior of the dipole-dipole scattering amplitude. 


• We found out that the unitarity bound for dipole-dipole cross section for very high 
energies is < 7 ( 7 * — 7 *) < ^ ln(l/a:). This result can be translated in the unitarity 

bound for 7 * — 7 * cross section after integration over ri^t and r 2 ,t in Eq. (HH). 

For Qi < Q? < Ql,at = we obtain: 


rrT,T{l* — 7*) 
(Tt,l{i* - 7*) 
(TlAi* - 7*) 


(TlAi* - 7*) 


< 


< 


< 


< 




ln(l/a;); 
ln(l/x); 
ln(l/a;); 



ln(l/a;); 


(6.43) 


In Eq. (I6.43|l for transverse polarized photon we used the logarithmic approximation in 
the integral over r^. Indeed, iT-pp oc 1/r/ and it should be integrated from 4/Q^^j to 
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4/Q^. Taking into acconnt that oc xG oc \/?hi.t>^ 2 ,t) in(i/*) (j;t. 21 j) i one can 

see that 


o'T,r(7* ■ 

-7 ) 

< 

Ct,t 

( 2 m ^)2 

\v?{l/x ); 

(6.44) 

- 

-7*) 

< 

Gt,l 

( 2 m ^)2 

3 

ln2(l/a;) ; 

(6.45) 

c^l,l(7* - 

-7*) 

< 

Cl,l 

( 2 m ^)2 

\vf{l/x ); 

(6.46) 


Therefore, only aT,T has the same energy dependence of the nnitarity bound as hadron 
-hadron cross section. 

• Our approach shows that the non-perturbative corrections at large 6 * should be taken 
into account in the Born cross section. Another way to treat this result is to say that 
the non-perturbative corrections can be taken into account only in the initial conditions 
as was discussed in Refs. | 2 HlEn|. We do not see that such corrections are needed in the 
kernel of the non-linear evolution equation [SEj as was argued in Refs. 12111221123 EH- 

The estimates at what energies such corrections will enter the game are presented and 
discussed. 

• Numerical calculations are performed for the value of the total cross section for acces¬ 
sible energies and virtualities. These predictions will be checked soon with new coming 
data. 

We hope that this paper will stimulate further experimental study of 7 * — 7* - processes 
which can give a very conclusive information on the saturation kinematic region in QCD. 
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